A Two-field Dilaton Model of Dark Energy 
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We investigate the cosmological evolution of a two-field model of dark energy where one is a 
dilaton field with canonical kinetic energy and the other is a phantom field with a negative kinetic 
energy term. A phase-plane analysis shows that the phantom-dominated scaling solution is the 
stable late-time attractor of this type of models. We find that during the evolution of the universe, 
■ the equation of state w changes from w > — 1 to w < — 1, which is consistent with the recent 

' observations. 

>• ' PACS numbers: 95.36.+X, 98.80.-k, 98.80.Es 

O ' 



O 

u 



o 

OS 

o 



- 1—^ 

X 



I. INTRODUCTION 



In recent years, observations of Type la supernovae (SNe la) [1, 2], cosmic microwave background (CMB) fluctua- 
tions [3, 4], and large-scale structures (LSS) [5, 6] indicate that the Universe is accelerating, therefore some form of 
• dark energy whose fractional energy density is about Qde = 0.70 must exist in the Universe to drive this acceleration. 
Dark energy has been one of the most active fields in modern cosmology since the discovery of accelerated expansion 
of our universe. Investigation on the nature of dark energy becomes one of the most important tasks for modern 
physics and modern astrophysics. Up to now, many candidates of dark energy have been proposed to fit various 
observations which include the simplest one, the Einstein's cosmological constant [7], or a dynamical scalar field, such 
as quintessence [8], phantom [9], k-essence [10], tachyon [11] and so on. The present data seem to slightly favor an 
evolving dark energy with the equation-of-state parameter (EoS) w < — 1 around present epoch and w > — 1 in the 
£NJ ■ near past. Obviously, w cannot cross —1 for quintessence or phantom alone. Some efforts have been made to build 
dark energy model whose EoS can cross the phantom divide. In a universe filled with quintessence and phantom 
fields this case can be realized easily. This implement of dark energy, called as quintom, has been first proposed in 
Ref. [12], where the quintom model with an exponential potential and the existence, stability of cosmological scaling 
solutions in the context of spatially homogeneous cosmological models have been investigated. Phase-plane analysis 
• ■ of the spatially flat FRW models shows that the phantom-dominated scaling solution is the unique late-time attractor 
and there exists a transition from w > — 1 to w < — 1 [13]. Wei and Cai [14] suggested a hessence model, in which 
a non-canonical complex scalar field plays the role of dark energy. The cosmological evolution of the hessence dark 
energy is also investigated; it is found that the big rip never appears in the hessence model even in the most general 
case while beyond particular potentials and interaction forms. 

The action of dilaton field in the presence of Einstein's cosmological constant has been derived in Ref. [15]. The 
potential is the counterpart of the Einstein's cosmological constant in the dilaton gravity theory. Since it can be 
reduced to the Einstein cosmological constant when the dilaton field is set to zero, the dilaton potential is called 
the cosmological constant term in the dilaton gravity theory. Compared to the ordinary scalar field, the action for 
phantom scalar field has only a sign difference before the kinetic term. Later, the explicit expression of the phantom 
potential have been given in Ref. [16]. A model of the Universe dominated by the dilaton field with a Liouville type 
potential has been presented in Ref. [17]. 

In this Letter, we investigate the cosmological evolution of a two-field model of dark energy where one is a dilaton 
field with canonical kinetic energy and the other is a phantom field with a negative kinetic energy term with Liouville 
type potentials. 
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II. EQUATIONS OF MOTION FOR THE TWO-FIELD DILATON MODEL 

Let us start from a 4-dimensional theory in which gravity is coupled to dilaton and Maxwell field with an action: 

S x = ( d 4 x^ (R - 28^(f) - Vi (</») + e-^F 2 ) , (1) 



where R is the scalar curvature, F 2 = F^F^ is the usual Maxwell contribution, a is an arbitrary constant governing 
the strength of the coupling between the dilaton and the Maxwell field, V\{<f>) is a potential of dilaton (f> which is given 
byRef. [15] 

W) = 0/1 2A 2 > 2 (3« 2 - IK** + (3 - a 2 )e 2a * + 8a 2 e^% (2) 
3(1 + a z ) 

here A is the cosmological constant. One can verify that the potential reduces to the Einstein cosmological constant 
when a = or <f> = 0. Compared to the action of the ordinary scalar fields, the phantom field has one negative kinetic 
term. In order to obtain a real action of the Einstein-Maxwell field in the presence of the phantom, we can make 
substitutions in the action as follows Ref. [16] 

4> — ► iip, a — > i/3, (3) 
where i is the imaginary unit. Thus we get the action 

S 2 = [ d A xJ—g (R + 2d^d^ - V 2 (VO + e~ 2 ^F 2 ) , (4) 



and the potential for the phantom field 

W) = o n 2 \ 2 ^ 2 (3/3 2 + IK** + (/? 2 + 3)e 2W - 8/3 2 e^-^]. (5) 
3(1 - P l ) 

One can also verify that, when (J = or ip = the action reduces to the Einstein-Maxwell action and when F 2 = 
the action reduces to the Einstein-phantom action. 

We consider the action in a simple model which contains a normal scalar field <f> and a negative-kinetic scalar field ip, 
assuming that there is no direct coupling between the phantom field and the normal scalar field with such potentials, 

S = Jd 4 x^g(R - 2d^<t> + 2d^d»ip - Vi(<l>) - V 2 (ip) + C m ) , (6) 

where C m represents the Lagrangian density of matter fields. Considering a flat Universe which is described by the 
Friedmann-Robertson- Walker metric, the homogeneous fields <j> and rp can be described by a fluid with an effective 
energy density p and an effective pressure P given by 

p = + ±vi(<f>) + \v 2 {^), (7) 

P = 4?-^- ±Vi(4i)- / 2 W- (8) 

The corresponding equation of state (EoS) parameter is given by 

i W )_ iy aW 

W = -. : 7 7 . (9) 

Then the equations of motion for the fields and the Friedmann equation can be written as 

3H 2 = n 2 (j? -i: 2 + Ivitf) + ^W) + P-yj , (12) 

where p 1 is the density of fluid with a barotropic equation of state P 7 = (7 — l)p 7 with 7 a constant and < 7 < 2 
(7 = 4/3 for radiation and 7=1 for dust mater). The equation (|12[) is the Friedmann constraint equation. 
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III. PHASE SPACE ANALYSIS AND THE CRITICAL POINTS 



In this section, we investigate the two-field Dilaton model via the conventional phase space analysis. Similar as in 
Ref. [18], we define the following new dimensionless variables 



X *~V3H' V4, ~ y/6H iVM* ^~ [^M]2 



= «yP7 

the equations of motion (fT0|) -(fl"2" |) can be rewritten as the following system of equations: 



^ = 3x {x\ - x\ + - l) + \^4>v1 > ( 13 ) 

= Zy<p (%l ~xl + ^z 2 ^j - ^X+xtyj, , (14) 

d\<p 
dN 



-x^x^ (r^ - 1) , (15) 

3x0, (xl -x% + - l) - 7;X,pyl , (16) 



dy^h ( 2 2 7 2 \ 1 

— ^ = 3y^, (z - x CT + -z J - -A^atyjty , (17) 



dXff, 
~dN 



-x^x^(r^-i), (is) 



o / 2 2 7 2 7 



where iV is the logarithm of the scale factor (N = In a), and the Fridemann constraint equation (11) becomes 

xl + yl-x%+yl + z 2 = 1. (20) 

Different from the case of a single exponential potential, the parameters and T here are variables of (f> and if). 
Strictly speaking, the above system is not an autonomous system. Thus, if we want to discuss the phase plane, we 
need to find the constraints on the potential, or equivalently the conditions under which the potential may have the 
property we require in order that we can get some explicit results. 

Critical points correspond to fixed points where ^ = 0, = 0, ^ = 0, ^ = 0, = 0, ^ = 0, jfc = 0. 
Observing these equations, one can find that the physically meaningful critical points (x^c, y<j>, c, X$ >c , X^ tC , y^, c ^ A^, jC ) 
of the system are: (Note that we will restrict our discussion of the existence and stability of critical points in the 
expanding universes with H > 0). 

(i) . (A^, c ^ 0, A^, c ^ 0) could be fixed by T = 1, = 1; 

(ii) . (x^c = 0, y\ c + y\ c = 1, A 0jC = 0, x^ >c = 0, A v ,, c = 0); 

(iii) . (x^c = 0, y^ c = 0, X^ c = any, x^ >c = 0, y</, )C = 0, X^ tC = any). 

In the case of (i), when 1^ = 1, r^, = 1, then X^^ — 2 -^, A^ c = and the two-fields potentials are given by 

F!^) = A(e"^ + 3), (21) 
V 2 (i>) = 2Ae- 2 ^. (22) 

Thus we can determine the important parameters of the dilaton field and the phantom field: 

a = or V3; and (3 — 1. 

v3 

In this case, the equations (fT3"]) - (|19p have one two-dimensional hyperbola (Type A) embedded in four-dimensional 
phase-space corresponding to kinetic-dominated solutions, with EoS w = 1 and fractional energy density Que = 1; 
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Type x$ y$ Xij, yj, z w SIde Stability 
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2 2 


= 1 


2 2-1 

x <t, ~ x ij> — 1 








1 
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unstable 
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1 













-1 + 


2 

'.'i K - 
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unstable 
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^/3k 2 7 (2-7) 















3 
2k 2 


unstable 


D 















-1 - 


2 
3k 2 


1 


stable 



TABLE I: The properties of the critical points in a spatially flat FRW universe containing a phantom field and a normal scalar 
field in the case of (i). 



a fixed point (Type B) corresponding to a dilaton-dominated solution, with w = — 1 + 2/3/t 2 and f2rjE = 1; a fixed 
point (Type C) corresponding to a fluid-dilaton-dominated solution, with w = and Ode = 3/2k 2 ; and a fixed point 
(Type D) corresponding to a phantom-dominated solution, with w = —1 — 2/3k 2 and J7de = 1 (listed in Table 1). 

In order to study the stability of the critical points, using the Friedmann constraint equation (|20p we can reduce 
Eqs.(fT3 |) - (|19|) to four independent equations. Substituting linear perturbations x<j> — > x^, + 5x<f,, y<p ^ y<p + Sy^, 
x^p — > Xip + Sx a and y^ — > + <$?/,/, into the four independent equations, we obtain the equation of perturbations to 



the first-order: 



dx'^ = 3 (Sx 2 ^ - x^ + ^z 2 - lj 6x,p + X^y^Sy^ - Gx^x^Sx^, (23) 

W$ = ^2a; - ^A^ Sxcf, + (x 2 -x 2 b + ^z 2 - lj Sy^ - Qx^y^Sx^, (24) 

5x'^ = Qx^x^Sx^, - 3 (x 2 ^ - 3x 2 -I- ^z 2 - l) 8x$ - X^y^Sy^, (25) 

Sy'^ = Qx^x^Sy^, + 3y c f > (~2x <j> - -\^)6x^ - 3 (ix 2 ^ - x\ + ^z 2 - lj Sy^. (26) 



The linear perturbations of system (|23j) -(|26 |) about each fixed point gives four eigenvalues. The theory of stability 
requires that the real part of all eigenvalues should be negative. So we have: 
Type A (the kinetic-dominated solution): 

mi = 3, m-i = 0, m3 = 3(2 — 7), 7774 = 3(1 ± _ ), 



V3k 



indicating that this solution is always unstable. 
Type B (the dilaton-dominated solution): 



6 6 „ 12 „ 

7771 = -T, 7772 = 7773 = — - 37, 7774 = —~ ~ 3, 
K z K z K z 



indicating that this solution is also unstable. 
Type C (the fluid-dilaton-dominated solution): 



777 - 37 777 - 37 3 777 - 3 ( 2 ~ ^ f| + L M^^) 

mi- Y , r77 2 _ Y - 3, 777 3 , 4 _ _M±yi_____ 



indicating that this solution is still unstable. 
Type D (the phantom-dominated solution): 

1 1 2 „ 

7771 = n , "72 = "13 = n - 3, 777 4 = <T - 37, 

K z K z K z 

indicating that this solution is stable. 

In the case of (ii) (x^,. c = 0, y^ c + y^ c = 1, A^ )C = 0, x^ )C = 0, At/, )C = 0), the equations (fTB|) - P^)l have one fixed 
point Type E embedded in six-dimensional phase-space which corresponds to eigenvalues (—3, —37, 0, 0, 0, 0), w = — 1 
and f^DE = 1- This indicates that the critical point is a de Sitter attractor. 

In the case of (hi) {x^ tC — 0, y$ tC = 0, X^ tC = any, x^ tC = 0, y^, tC — 0, A^ ;C = any), the equations (|13p - (fT9|) 
have one fixed point Type F embedded in six-dimensional phase-space which corresponds to eigenvalues (—3(7 — 
2)/2, —37/2, 0, —3(7 — 2) /2, —37/2, 0), w becomes meaningless and Ode = 0. This indicates that the critical point is 
not a dynamical attractor. 
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FIG. 1: The evolution of the equation of state w of the two-field dilaton model of dark energy. The blue line represents the 
EoS of the two-field model, the black line represents the single dilaton model and the red line represents the single phantom 
model. 



IV. NUMERICAL STUDIES 

Our numerical studies indicate that the EoS parameter w changes from > —1 to < —1 as shown in Figure 1. We 
have assumed that there is no direct coupling between the phantom field and the normal scalar field in this paper. 
Without the loss of generality, the initial conditions <fi(0), ip(0), <p(0) and V>(0) can be fixed in order to get the EoS 
today (w — —1.02, a — 1) [19], and the energy density of dark energy today Qde = 0.70. The blue line represents 
the EoS of the two-field dilaton model, the black line represents the single dilaton model and the red line represents 
the single phantom model. 

V. CONCLUSIONS AND DISCUSSIONS 

In summary, we have investigated the possibility of constructing a two-field dark energy model which has the 
equation of state w crossing —1 by using the dilaton and phantom fields. We have made a phase-space analysis of the 
evolution for a spatially flat FRW universe filled with a barotropic fluid and phantom-dilaton fields. It is shown that 
there exists the stable late-time attractor solution in the model. Also, we showed that the equation of state w can 
cross —1 naturally. So the two-field dilaton field is a viable candidate for dark energy. 

It is apparent that our model is also plagued with the instability problem at the quantum level which makes its 
existence doubtful. In fact, this is a common problem for nearly all phantom models. However, as argued by Carroll 
et al. [20], these models might be phenomenologically viable if considered as effective field theories valid only up to 
a certain momentum cutoff. According to their discussions, the instability timescale of the phantom quanta can be 
greater than the age of the universe provided that the cutoff is at or below 100 MeV. In this sense, the phantom 
quanta are stable against decay into gravitons and other particles. Therefore, considering astronomical observations 
favoring the phantom model for dark energy, it remains open if the phantom matter exists and acts as dark energy. 
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